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On the Origin of Broad Fe Ka & Hi Ra Lines in AGNfj 
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ABSTRACT 

We examine the properties of the Fe emission lines that arise near 6.4 keV 
in the ASCA spectra of AGN. Our emphasis is on the Seyfert 1 galaxies where 
O ! broad and apparently complex Fe Ka emission is observed. We consider various 

origins for the line but focus on the pros and cons for line emitting accretion disk 






OS 



(N 



models. We develop a simple model of an illuminated disk capable of producing 
both X-ray and optical lines from a disk. The model is able to reproduce the 
observed Fe Ka FWHM ratio as well as the radii of maximum emissivity 
implied by the profile redshifts. The overall profile shapes however do not fit 
well the predictions of our disk illumination model nor do we derive always 
consistent disk inclinations for the two lines. We conclude that the evidence for 
and against an accretion disk origin for the Fe Ka emission is equal at best. 
The bulk of the data requires a very disparate set of line fits which shed little 
light on a coherent physical model. We briefly consider alternatives to disk 
emission models and show that a simple bicone model can reproduce the Fe 
line profiles equally well. 
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1. Introduction 



The idea that a significant fraction of the broad line emission in Active Galactic 
Nuclei (AGN) comes directly from an accretion disk is attractive for several reasons beyond 
the obvious one of proving the black hole paradigm. The disk provides a high density 
medium where the integrated intensities of the optical Feii lines can be accounted for. 
An illuminated disk can also more easily produce the strong Balmer lines (i.e. the La/H/3 
problem: see Collin- Souffrin & Dumont 1990). Apparent density and kinematic differences 
between high and low ionization lines (Gaskell 1982; Sulentic 1989) can also be explained 
if the latter arise from a disk. Naturally there was considerable excitement when the 
double-peaked Balmer lines in Arpl02B were discovered (Chen, Halpern and Filippenko 
1989). The line profiles in that AGN could be well-fit by the predictions of a model for line 
emission from a relativistic Keplerian disk (Chen and Halpern 1989). The fit to Arpl02B 
was quite good but the statistical implications were not. The best disk model fits to lines 
in Arpl02B imply an intermediate viewing angle and that leads to the expectation of many 
sources with double-peaked emission lines. A search for more double peaked and peculiar 
profiles among radio-loud AGN revealed only a handful that could be reasonably well fit by 
disk models (Eracleous and Halpern 1994). 

AGN with optical line profiles like Arpl02B are apparently quite rare. They also 
appear to be almost uniquely radio-loud sources. Either most radio-loud disk emitters (and 
all radio-quiet ones) produce emission at larger radii where double peaked lines are not 
expected or the disk contribution to the broad line emission is small or negligible. The 
unsatisfactory alternative would be that we are viewing virtually all AGN aligned with the 
disk axis near pole-on. A comparison between line profile shifts/asymmetries with disk 
model predictions (Sulentic et al. 1990) suggests that at least part of the line emission 
in most quasars must arise from another source. This conclusion is based on the large 
number of line profiles that show blue shifts and/or blue asymmetries (implying that the 
red side of the profile is brighter) in conflict with disk model predictions. Even the most 
promising double-peaked sources (e.g. Arpl02B, 3C390.3) require fine tuning of the disk 
models (see e.g. Gaskell 1996; Zheng et al. 1991) because the profile peaks appear to vary 
out of phase. A new double peaked source arose in Pictor A with first a red, and then a 
blue, peak developing over a five year period (Sulentic et al. 1995a). Such variations and 
the overall rarity of double-peaked profiles appear to be more easily explained with models 
invoking a predominance of radial motion. Such models also have the advantage that they 
can directly account for the rarity of peculiar line profiles because such profiles would arise 
at an extremum in the orientation/kinematic domain (Sulentic et al. 1995a). 

A possible kinematic difference between high and low ionization lines (HIL and LIL; 
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HIL include CivA1549 and Heii lines; LIL are Balmer lines, FEH opt lines, and MgiiA2800) 
was one of the reasons that models for a two zone broad line region (disk + something 
else) became necessary and popular. The evidence revolves around a possible systematic 
blueshift of the HIL with respect to the LIL (the LIL showing centroid velocities similar 
to the source rest frame). The most peculiar line profiles were found in the radio-loud 
AGN population and subsequent searches for disk candidates focused there (Eracleous and 
Halpern 1994). Ironically the largest sample yet studied for line profile shifts indicates that 
it is only the radio-quiet AGN that show a systematic HIL blueshift with respect to the LIL 
(Sulentic et al. 1995b; Marziani et al. 1996). 

Recent advances in X-ray spectral resolution offered by the Advanced Satellite for 
Cosmology and Astrophysics (ASCA; Tanaka, Holt, & Inoue 1994) provide a possible 
new form of evidence for disk emission in (predominantly) radio-quiet AGN. We may be 
observing the very innermost regions of the disk that could be producing an iron Ka line 
by Compton reflection or photo-ionization (e. g. Tanaka et al. 1995; Fabian et al. 1995; 
Dabrowski et al. 1997; Fanton et al. 1997). Strong and broad line emission is often observed 
near 6.4 keV in AGN spectra (Mushotzky, Done, & Pounds 1993; Mushotzky 1997; Nandra 
et al. 1997a,b, and references therein). Figure 1 illustrates the difficulties that we face in 
quantifying the X-ray line emission. While ASCA provides a tremendous advance in X-ray 
spectroscopy, the resolution (2-5% at 6 keV) and S/N obtainable still lag far behind optical 
data. We show a typical modern observation of the Ha line profile for Arpl02B degraded 
to the resolution and S/N of a typical ASCA spectrum. Narrow lines near Ha arising 
from [Nil] AA6548, 6583, [OI]A6300 and [SII]AA6717,6730 were removed from the spectrum 
before blurring because the Fe Kq emission in AGN has no known analogous contaminants, 
albeit Fe Ka might involve a blend of multiple, broad and narrow, components. While 
optical observations can discriminate between a disk model and a broad gaussian, ASCA 
observations cannot (see Fig. 1). 

We consider Fe line detection statistics for different AGN classes in § |2|. In § 3 we 
consider an accretion disk origin for the Fe Ka and Hi Ha lines. In § 3 we model an 
illuminated disk (e.g. Dumont & Collin-Souffrin 1986) to account for both optical and 
X-ray emission lines simultaneously. In § 4 we summarize problems with the accretion disk 
interpretation for Fe Ka and consider an alternative model for the X-ray line emission. 



2. X-ray and Optical Observations 

The energy of the emitted Fe Ka depends on ionization stage. The ASCA resolution 
allows us to distinguish line emission from three main ranges of ionization: (1) emission 
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from neutral or weekly ionized iron ("cold iron" emission; ionization stage <FEXVIIl), as 
expected from cold, dense gas in an accretion disk; (2) emission from mildly to highly 
ionized iron at 6.45 keV £ energy ,$ 6.7 keV, for Fexix< ionization stage <Fexxv; 
(3) emission from H-like Fe at 6.97 keV (Fexxvi). The most robust measurement is the 
centroid energy of the line, but even this can become confused when the line is broad. 
Determination of the profile width, asymmetry index or equivalent width are more difficult 
and model dependent. These parameters, and the derived y 2 goodness-of-fit are quite 
sensitive to the continuum fit which includes a basic power law modified by both emission 
and absorption processes (correlated to the line itself). 



2.1. Statistics of Emission lines near 6.4 keV 

Table 1 summarizes the published Fe line parameters for most AGN with detected 
Fe Ka emission in ASCA spectra. The list should be reasonably complete through 1996. 
The format for Table 1 is as follows: Column 1 - abbreviated (2000) position taken from 
the 7th edition of the Veron & Veron catalog; Column 2 - a common name for the source 
as used in the ASCA-related literature; Column 3 - AGN type; Column 4 - optically 
determined redshift; Columns 5, 6, and 7 - the derived Fe Ka line peak centroid energy, 
profile Gaussian width (keV) and EW (eV), respectively; Column 8 - reference(s) to the 
source of the Fe Ka parameters given in the three previous columns. Table 1 reflects the 
best currently derived line parameters. In some cases a range is given to reflect a range 
of reasonable solutions, providing additional insight into the uncertainties associated with 
these numbers. References should be consulted for details of adopted continuum and line 
properties. We focus our discussion on the Seyfert galaxies because: (1) that is where the 
largest body of X-ray emission line data exists and the statistics of Fe Ka emission from 
broad-line radio galaxies and quasars are not yet clear and (2) the Fe Ka spectra of Sy 1 
galaxies sometimes show "disk-like" profiles. 

A recent study by Nandra et al. (1997a,b) has collected much of the data and reduced 
it uniformly. The best model fits to Seyfert 1 spectra always show Fe emission with 
approximately 80% of the detections better fit by a resolved broad line. Nandra et al. 
(1997b) find a distribution of profile (Gaussian fit a) widths with: (1) about half of the 
sample clustered at 0.7±0.1 keV and (2) the other half of the sample extending continuously 
from 0.5 keV down to unresolved instrumental width (~0.1 keV). The profiles frequently 
show a narrow component centered very near to 6.4 keV accompanied by an asymmetry 
towards the red. In several cases the spectra suggest a multicomponent structure within 
this red wing rather than a smooth extension. In most cases it is impossible to determine 
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if this multicomponent structure is real or due to instrumental effects and low S/N. The 
distribution of line equivalent widths extends from 50-600 eV with a mean EW~160. 
Typical uncertainties for the individual estimates cover most of this range. The major point 
of consensus is that we typically observe a broad component extending 1-1.5 keV to the red 
of a narrow peak that is centered very close to 6.4 keV. The two simplest models of the 
line involve: (1) a single asymmetric feature (§ 3) or (2) a combination of two symmetric 
features one narrow/unshifted and another broad/redshifted (§ 4). 

The situation for radio quiet quasars is very confused with many showing no evidence 
of rest frame emission near 6.4 keV. At the same time sources from 3C273 to S5 0014+813 
(Elvis et al. 1994) at z=3.384 show evidence for a line. Rest frame centroid energies for 
these detections range from 6.4-6.9keV. Few observations for radio-loud sources have been 
published and the situation is not clearly defined. No radio-loud sources as yet show 
broad Fe profiles like IRAS18325-5926 or Sy l's, as further discussed in § |]|. In 3C390.3 
(Eracleous et al. 1996) the Ha and Ka widths are similar while for 3C120 Fe Ka may 
show either a single broad line or 3 individual narrower features. A recent observation for 
Pictor A (Halpern et al. 1997) shows no Fe line. Recent claims (Reynolds 1997) that radio 
loud sources show much broader Fe Ka profiles than radio quiets are based upon very 
broad and uncertain detections for 3C120 and 3C382. 

Seyfert 2 sources usually show narrow Fe Ka emission. The results from a large 
GINGA sample (Smith & Done 1996) indicate that most lines have a centroid near 6.4 keV 
with a subset centered between 6.6-7.0 keV. In two of the best studied cases (Mark 3 and 
NGC1068) there is evidence for both neutral and ionized Fe lines, possibly components of 
Feu emission from all three ionization ranges listed above. The most unusual Sy2 sources 
are IRAS 18325-5926 and MCG-5-23-16 where very broad lines were detected. The former 
source shows a broad (spanning 4.5-7 keV) profile, red (low energy) asymmetry and narrow 
peak on the blue (high energy) side centered at 6.8 keV. It is unclear whether this is a single 
line or a blend of several features but it is reasonably certain that the extended emission 
redward of 6.4keV is real. MCG-5-23-16 shows a broad component extending from 5-8 
keV with a narrow peak near 6.4 keV. Multiple components are again a distinct possibility. 
The "hidden" quasar IRAS 09104+4109 (Fabian et al. 1994b) has been proposed as an 
ultraluminous example of the Sy2 class. It shows a very strong but narrow Fe line centered 
at 6.65 keV. 
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2.2. BLR Ha Statistics and Comparison with Fe Ka 

We have obtained, or gained access to, optical spectra for almost every Seyfert 1 
galaxy that shows a strong Fe Ka (ASCA) detection so far. Observations that have not 
been presented elsewhere are summarized in Table 2. The format is as follows: Column 1 
- abbreviated IAU designation; Column 2 - other name; Column 3 - date of observation; 
Column 4 - universal time of observation; Column 5 - exposure time in seconds; Column 
6 - observatory; Column 7 - telescope; Column 8 - spectrograph; Column 9 - grating and 
Column 10 - resolution FWHM (in A). 

Table 3 summarizes the Hi Ha profile properties for both old and new observations. 
The format is as follows: Column 1 - IAU code name; Column 2 - source name; Column 3 - 
Seyfert type; Column 4 - EW in A; Column 5 - FWHM in km/s; Column 6 - uncertainty at 
2xo" confidence level; Column 7 - centroid line shift (at half maximum) in km s" 1 ; Column 
8 - line centroid uncertainty; Column 9 - reference code. Uncertainty of Ha EW values 
should be ~ 10%. Figure 2 shows a comparison of Ha and Fe Ka profile FWHM values. 
The values plotted for Fe Ka come from the (2.35 x) a values given in Table 1. They were 
derived from the broad line fits to the spectra (see Table 1 references). They were usually 
estimated by fitting the X-ray continuum with (a) a nonthermal power-law, (b) a soft X-ray 
absorption law and (c) a hard reflection component. In general inclusion of component (c) 
enhances the continuum underlying the line resulting in a reduction of the derived EW and 
FWHM values. In some cases plausible solutions can cause the line to become unresolved 
(e.g. MCG-02-58-22). 

Figure 2 shows that Fe Ka line is systematically broader than Ha in all 
radio-quiet AGN. In the case of MCG-6-30-15 (dots in Figure 6) and NGC 4151, 
FWHM(Ka)/FWHM(Ha) « 40 and 10 respectively. The limited data for radio-loud 
sources shows a smaller difference in FWHM. For the Seyfert 1 galaxies listed in Table 3 we 
find: (1) the hypothesis that Fe Ka and Ha FWHM come from the same distribution is 
rejected at a confidence level £ 0.995 using generalized (including upper limits for Fe Ka) 
Wilcoxon tests and (2) Fe Ka and Ha FWHM are not significantly correlated. A Kendall 
test with the inclusion of censored data yields a probability P~ 60% that a correlation 
is absent. EW values of Fe Ka and Hi Ha are also not significantly correlated. The 
Spearman correlation coefficient is rg « 0.31. For 25 data points the probability to have this 
value in a random sample of observations taken from an uncorrelated parent populations is 
12 %. Restricting attention to the objects studied by Nandra et al. (1997), we do not find 
any significant correlation for half maximum centroid shifts. The most interesting trend 
emerging from a lineshift comparison is that the objects with a strongly redshifted broad 
component in Fe Ka show symmetric (unshifted) and relatively narrow Hi Ha line profiles. 
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2.3. Possible Sources of Fe Ka Line Emission 

Our considerations of origins for Fe Ka are focused on the Seyfert galaxies where 
the bulk of the Fe Ka spectral detections of reasonable S/N are found. In a unification 
driven approach we assume that: (1) a supermassive black hole and associated accretion 
disk lie at the center of each AGN, (2) Seyfert 2 galaxies (at least some of them) are Seyfert 
l's seen edge-on (i £ 45°) and (3) the optical broad line emission region in Seyfert 2's is 
obscured by the putative cold dust torus. In this framework we can identify at least six 
possible origins for Fe Ka line emission (listed roughly in order of decreasing distance from 
the black hole): (A) cold Fe Ka emission from the obscuring dust torus which is assumed 
to be optically thin at ~ 6keV in most cases (see e.g. Mushotzky et al. 1993: Ghisellini, 
Haardt & Matt 1994; Weaver et al. 1996); (B) Fe Ka emission from the classical BLR; 
(C) Fe Ka emission from a photoionized scattering region (a warm absorber); (D) cold 
Fe Ka emission due to fluorescence reflection (or emission) from the accretion disk; (E) 
Fe Ka emission from highly ionized coronal gas. Coronal gas (T £ 10 9 K) can give rise to 
strong Fe Ka emission at 6.67 or 6.97 KeV; (F) cold Fe Ka emission from a (non-disk) 
cloud distribution (not necessarily responsible for optical/UV LIL and/or HIL). Option 
(A), (B), (C), and (E) are briefly discussed below, as they appear rather unlikely to be the 
dominant source of Fe Ka emission. Option (D) and (F) deserve a much deeper discussion, 
deferred to § and § 4 respectively. Our model considerations focus on MCG-6-30-15 as 
the best studied source with a characteristic broad Fe Ka line profile. The broadband 
spectrum of MGC-6-30-15 is consistent with a typical AGN continuum as described by 
Mathews & Ferland (1987; MF87) except that the sub-millimeter break may occur at higher 
energy than implied by their continuum ( £ lxlO -3 Ryd in MCG-6-30-15). We nonetheless 
assume that the sub-millimeter break occurrs at lxlO -5 Ryd as parameterized by MF87. 
The ionizing luminosity of MCG -6-30-15 is Li on « 6 x 10 43 ergs s _1 from logz/L(i/) « 29.0 
ergs s _1 at v « 14 Hz (Ho = 75km s _1 Mpc -1 , and qo = is assumed throught the paper; 
the corrected heliocentric radial velocity of MCG-06-30-15 measured on our spectrum is 
v r h « 2317 km s _1 , which implies a distance ~ 30.9 Mpc.) 

Table 4 reports expected line luminosities of MGC -6-30-15 for cases A, B. C, E. 
In Column (2) we report the Fe Ka luminosity computed in the case of gas with solar 
chemical abundances, while in Column (3) we report the Fe Ka luminosity in the case of a 
factor 10 iron overabundance with respect to the solar value. In Column (4) the dominant 
ionization stage of iron is reported. Column (5) lists the expected Hi Ha luminosity. 
Luminosities are in ergs s _1 . These values should be compared to an observed total Fe Ka 
luminosity logL(FE Ka) ~ 41.2, and to logL(Hl Ha) « 40.8 . The L(Ha) value has been 
corrected for galactic extinction assuming Ab ~ 0.15, but not for internal absorption. The 
Fe Ka narrow component luminosity is ~ 0.3-0.5 of the total L(Fe Ka). 



2.3.1. Emission from a Dusty Torus 

Option A would only give rise to a narrow Fe Ka line centered near rest energy. Strong 
variations observed in the "narrow" 6.4 keV part of the MCG-6-30-15 spectrum complicate 
this interpretation (Iwasawa et al. 1996b). For MCG-6-30-15 there is also a photon deficit 
problem. Table 4 reports the results of explorative CLOUDY (Ferland 1996) photoionization 
computations. A simple, optically thin model of a dusty torus, located between 3 and 10 
pc, and photoionized by the central continuum source, produces a Fe Ka luminosity that 
is more than a factor ~ 100 below the observed Fe Ka luminosity. This remains true even 
if the attempt to ascribe to torus emission the "narrow component" only of the Fe Ka 
line. The situation may be markedly different in Sy 2 galaxies, where a narrow Fe Ka line 
is observed against a strongly absorbed X-ray continuum (Weaver et al. 1996; Ghisellini, 
Haardt, & Matt, 1994; Krolik, Madau & Zycki 1994). 



2.3.2. Warm Absorber 

The warm absorber (option B) is also unlikely to be the dominant source of Fe Ka 
line emission (Reynolds & Fabian 1995; see Table 4). The photoionization calculations 
presented here are based on the model by Reynolds & Fabian: gas with Tsa 10 5 K, n e ~ 10 6 
cm -3 , and moderate column density N c ~ 10 22 cm -2 , located within the BLR or at the outer 
edge of the BLR is photoionized by the AGN continuum. This model is not the only one 
proposed for the warm absorber (e. g., Mathur, et al. 1994), but it is probably adequate 
for low luminosity radio-quiet AGN. There are two main obstacles to "warm iron" Fe Ka 
emission, (a) a photon deficit; (b) the effects of Compton scattering on the line profile 
(Fabian et al. 1995). As a central corona is likely to be Compton thick, extreme broadening 
may occur (~ 1 keV; Fabian et al. 1995; Nandra et al. 1996). 



2.3.3. Classical BLR 

The main problem with Fe Ka BLR emission are that (1) there is a flux deficit, unless 
high iron abundances are invoked (see Table 4), and (2) Fe Ka and Balmer line width 
are inconsistent at least in Seyfert 1 galaxies (§ |2.2| ). Both problems are solvable if only 
the narrow component of the Fe Ka line is ascribed to the BLR. Our photoionization 
simulations show that, if iron abundance is 10 times the solar value, then the luminosity 
of Fe Ka produced within the BLR agrees with the observed luminosity of the Fe Ka 
narrow component. 
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The so called "X-ray Baldwin effect" (anticorrelation between Fe Ka EW and 
X-ray luminosity) has been interpreted as evidence for a common origin of Fe Ka and 
CivA1549 emission in the BLR (Iwasawa & Taniguchi 1993). The ASCA data do not 
support such an anticorrelation. The EW of Ka is strongly dependent on the fit of the 
X-ray continuum. If we consider the Iwasawa & Taniguchi (1993) objects for which ASCA 
observations are available (16 AGN, mainly low-luminosity), and that were fitted with a 
reflection component whenever appropriate, we find that the anticorrelation disappears (a 
Kendall's r, generalized for the inclusion of censored data, yields ~ 15 % probability of no 
correlation). We must point out, however, that if we substitute the ASCA data, and keep 
the other Iwasawa & Taniguchi (1993) values (a total of 38 objects), obtained from GINGA 
observations, then the correlation does not go away (Kendall's r =>- P m 0.003-0.007). As 
the resolution of the GINGA data was substantially lower, and the typical profile in Sy 1 
galaxies resembles the one observed in MCG -06-30-15, with a broad redshifted component 
and a narrower component at 6.4 keV, a significant part of the broad companent may have 
been lost in the fit of the underlying continuum. The apparent lack of any "X-Ray Baldwin 
effect" in the ASCA data is therefore consistent with the idea that the Fe Ka narrow 
component only is emitted within the BLR (Marziani, et al. in preparation). 



2.3.4. Hot Corona 

A contribution from "hot iron" (more than 17 times ionized) cannot be ruled out for 
Seyfert nuclei because a high energy residual is observed for most Fe Ka disk candidates 
(see discussion in Section 4). We will show that, in fact, hot iron emission is expected from 
a photoionized disk. However, the most robust result of ASCA observations is that hot 
gas cannot be the dominant source of Fe Ka emission in Sy 1 galaxies, as the line energy 
is almost always around 6.4 keV, implying low ionization for iron (Nandra et al. 1997a). 
Nevertheless, coronal gas at a temperature close to the virial temperature T vir ^5x 10 11 
M 8 rjfi K, within r £ 10 -100R g (where the gravitational radius R g = GM/c 2 ) should be 
a strong Fe Ka emitter. The resulting Fe Ka luminosity is strongly dependent on the 
amount of matter that makes up such a corona. In Table 4, we compute a model assuming 
T=10 9 K, and N c fa 10 25 cm -2 . Redshifted hot coronal gas in addition to BLR emission 
provides an alternative to accretion disk emission. The observed profile may be reproduced 
equally well. The combined effects of gravitational redshift and increasing rest energy can 
produce centroid values }t 5.0 at 6 R g (6.19 keV at 20 R g ) if the ionization stage of iron is 
high. However, this requires some sort of fine tuning. 
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3. Broad Fe Ka and HI Ha Accretion Disk Models 

Two independent lines of recent work consider emission from an accretion disk for 
Fe Ka (Zycki & Czerny 1994; Matt, Perola & Stella 1992; Fabian et al. 1995; Dabrowski et 
al. 1997) and for Hi Ha (Dumont & Collin-Souffrin 1990a,b & Collin- Souffrin & Dumont 
1990; Rokaki, Collin-Souffrin & Magnan 1993, Eracleous & Halpern 1994). The first order 
fits to the Fe Ka line are qualitatively more satisfactory than fits to optical profiles. Model 
profiles computed for a rotating disk surrounding a Schwarzschild or Kerr black hole are 
most clearly consistent with the Fe Ka line profile observed in MCG-6-30-15 (dashed 
profile in Figure 6 shows Tanaka et al (1995) fit). However, as the number of Fe Ka 
detections increase, and line variability measures become more accurate, the disquieting 
trend is that model fits become poorer and require an ever larger range in the parameter 
space for disk models (e.g. Iwasawa et al. 1996b; Weaver et al. 1997). There is no evidence 
for convergence toward some characteristic disk model. Only a few other sources (Fairall 9, 
NGC4151, NGC5548, IC4329A) show structure unambiguously similar to MCG-6-30-15. 
This may be due to real profile differences or to the large variations in S/N of the data. 
The latter possibility is likely because the composite profile constructed by Nandra et al. 
(1997) from about 16 objects (see solid curve in Figure 6) is qualitatively similar to the 
MCG-6-30-15 profile. To summarize, one can fit almost any Fe Ka line profile with a disk 
model especially: (a) when one fits low S/N X-ray line profiles to models where most of the 
emission is produced at the innermost disk radii and (b) when the entire parameter space 
of Schwarzschild and Kerr models is available. Coupling this with assumptions about an 
emissivity law allows one to fit a disk model to any profile. This implies that a greater 
degree of self consistency is needed. 

In the optical, where the lines are better defined the problems are basically the same, 
with the additional difficulty of very poor agreement between models and observations. 
Objects like Arpl02B are rare-forcing disk emission advocates to argue that (1) objects 
with similar double-peaked profiles are most often obscured, and/or (2) that most emitters 
produce the line at larger radii in the disk (making the discovery of Arpl02B somewhat 
miraculous). Truncation of line emission at r ^ 1000 R g demands a physical motivation 
that has yet to be found (see Eracleous & Halpern 1994). 

Nonetheless, an accretion disk origin for the Fe line emission might offer a 
straightforward explanation for some important observational findings. If Fe Ka is emitted 
extremely close to the central black hole then that line emission is most sensitive to disk 
structure. This might explain why a number of intermediate redshift (z ~ 1) quasars show 
neither strong Fe Ka nor a strong reflection continuum component (Nandra et al. 1997c). 
Disk structure could be markedly different in radio loud AGN as well. The standard a-disk 
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theory suggests that the region dominated by radiation pressure may not exist if the 
dimensionless accretion rate hi £ 0.4 (Shakura & Sunyaev 1973; Frank, King, & Raine 
1992). This case corresponds to the lowest luminosity AGN, which are thought to radiate 
in the sub-Eddington regime. In intermediate redshift quasars, the dimensionless accretion 
rate may approach unity, invalidating the thin disk assumption, and the central regions of 
the disk may puff up (Abramowicz, Calvani, & Nobili 1980; see also Calvani, Marziani, & 
Padovani 1989 for a review). In this case, we expect that Fe Ka is emitted only in the 
narrow funnel of the inflated disk, as its external layers would become too cool (Madau 
1988). This would explain the rare detections of Fe Ka in intermediate and high redshift 
quasars, including both the radio-loud and radio-quiet variety, and observed dependence on 
luminosity of the line profile (Nandra et al. 1996; Nandra et al. 1997c). If Fe Ka is emitted 
in the funnel, then it can only be seen if the disk is observed nearly face-on, the rarest of 
the occurrences. It is intriguing that the only quasar where strong Fe Ka is detected, PG 
1116+215 (Nandra et al. 1996), fits very well with this scenario. Model fitting of Fe Ka 
suggests i ~ 0°. The very narrow Hi Ha profile, and the extremely strong FEII opt emission 
make it a close replica of I Zw 1. These features can be explained in terms of emission from 
a face-on disk (Marziani et al. 1996). 

LIL, HIL and Fe Ka might all arise from the same component with X-ray, UV and 
optical line emissivities peaking at different radii from the central continuum source. The 
FWHM difference described in § 2.3 could be accounted for if the X-ray lines were produced 
in a disk with Fe Ka produced near the center and the Balmer lines further out. At first 
glance the lack of a correlation between Fe Ka and Ha FWHM appears to pose a problem 
for models of this kind. This is not necessarily the case however because in disk models the 
width of Ha (FWHM) is driven by illuminating flux while FWHM Fe Ka is sensitively 
driven by the inner disk emitting radius. 



3.1. A Disk Illumination Model for Both Ha and Fe Ka 

In MCG -06-30-15, the Eddington ratio is r] Edd = L bo i/L Ed d £ 0.34 (M»2x 10 23 g 
s _1 ), consistent with the assumption of a thin disk down to the innermost stable orbit. We 
estimate a lower limit to the mass of the central black hole Mbh ~4x 1O 6 M from the 
assumption of cloud virial motion, following Padovani (1989). We assume M=10 7 M Q and 
M pa 2 x 10 23 g s _1 in the computation of the disk structure (if line emission takes place in 
a disk the mass estimate should be divided by sin(< i >) pa 0.5). We assume that this is 
appropriate for MCG-6-30-15 and all low-luminosity Seyfert l's studied in this work. 



We model the region where the continuum is produced as a sphere of radius r pa 100 R. 



k 
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with r cs ps 0.25. Radiation emitted in this region directly illuminates the disk, and is also 
partly scattered towards the disk by the free electrons in the sphere. In addition, a spherical 
halo of free electrons may surround the disk, up to 10 6 R g , again with r es ~ 0.25. Given the 
large uncertainty in the ionizing luminosity of AGN, we compute models for two extreme 
choices of the continuum: (a) the observed continuum, according to MF87. This choice 
probably overestimates the ionizing luminosity, as it includes the big blue bump and the 
reflection component that are thought to be produced within the disk; (b) a broken power 
law with slope —1.4 between lxlO 15 Hz and 50 keV, and —2.5 for energy £ 50 keV, which 
may yield an underestimate of the ionizing photon flux. Computation of the scattered 
flux as a function of disk radius is as carried out by Rokaki et al. (1993). Relativistic 
effects (most notably returning radiation i. e., radiation emitted by an accretion disk which 
returns to its surface due to gravitational focusing or the shape of the disk; Cunningham 
1976) are neglected in the scattering process, even if they are thought to be relevant in the 
innermost region of the disk. Returning radiation may enhance Fe Ka emission by a factor 
ps 2 (Martocchia & Matt 1996; see also Dabrowski et al. 1997). 

We consider here the radial disk structure appropriate for a Kerr metric, with the 
angular momentum per unit mass a = 0.998 M (c = G = 1; as shown by Thorne 1974, 
the radiation emitted by the disk and swallowed by the hole prevents spin up beyond 
a limiting state a/Mii m ps 0.998 M; the corresponding innermost stable orbit is at rps 
1.25 R g ) and use the equations for the radial and vertical structure of a standard a-disk 
model, with the corrective terms given by Novikov and Thorne (1973). The disk is divided 
into concentric annuli of fixed fractional width Ar/r ps 0.2. Surface density, column 
density, and illuminating flux are assumed constant within each annulus. We approximate 
the hydrostatic equilibrium (leaving the vertical structure unperturbed by the incoming 
radiation) assuming that the density within the disk increases linearly with column density. 
This is a crude approximation since the vertical structure of the disk should be recomputed 
to take into account the heating provided by the illuminating flux. In the outer layers of the 
innermost regions, if the density is too low and the radiation density too high (for values 
of £ = L ion /r 2 n c ~ 10 4 ergs s _1 cm), the gas is heated to T ~ 1 x 10 10 K. In this case the 
disk layer actually becomes part of a corona, and we assume that radiation penetrates to 
a layer within the disk of optical depth larger by Ardisk = 1, with flux reduced by a factor 
exp(— Araisk)- We employed the photoionization code CLOUDY (Ferland 1996) to compute 
the emissivity in the Fe Ka and Balmer lines. 



3.2. Illumination Model Results 
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3.2.1. Fe Ka 



We identify four vertical zones within the disk, (a) Very hot gas - a halo that has 
undergone a thermal runaway due to heating by scattered radiation in the lowest density 
layers of the disk (r < 1). (b) Warm gas in an "intermediate" region of higher density, 
whose electron temperature varies from £ 10 7 down to a few 10 4 K. (c) Ionized gas at 
a few 10 4 K; (d) cold, almost neutral gas which is less heated by the incoming radiation 
(T £ few x 10 3 K). Reflection models usually suggest that the "cold" Fe Ka emission 
is produced in the two deepest layer. In Table 5 we report the results of our calculations. 
In Column (1) we give an identification number for each model. In Columns (2) (3), (4), 
(5) the inner and outer radii, and the electron scattering optical depth for the inner sphere 
and the outer halo are reported. In Column (6) we list the shape of the continuum (as 
described in the previous section). Column (7), (8), and (9) yield the logarithms of the 
Fe Ka luminosity from cold iron emission, the ratio L(Fe Ka)n t/L(FE Ka)c id, and the 
log(L(Hl Ha)) respectively. 

Our model of the vertical disk structure allows us only a coarse estimate of the 
relative flux contributions from "cold" and "hot" iron. The computed ratio FEnot/FEcoid 
(s=s 0.1 — 0.3) is sensitive to returning radiation which is not taken into account. The 
predictions of our model for FEH O t/F E c id are nevertheless in qualitative agreements with 
those of Matt, Fabian & Ross (1996). We are also aware that mean escape probability 
formalisms are not adequate for the treatment of Compton scattering and radiative transfer 
of the Fe Ka line. Nevertheless, our photoionization calculations for cold Fe Ka agree 
with the results of Monte Carlo simulations of Fe Ka emission from a homogeneous layer 
of cold matter by George & Fabian (1991). They find that ~ 6% of the photons between 9 
and 30 keV are converted into Fe Ka photons. This yields L(Ka) ~ 1.2 x 10 41 ergs s _1 for 
MCG-6-30-15 (the observed total Fe Ka luminosity is L(Fe Ka) « 1.6 x 10 41 ergs s" 1 ). 

From Table 5 we deduce also that the observed luminosity values for MCG -6-30-15 
are in the range predicted by our disk models. We must remark that the choices of the 
continuum are probably extreme (new observations suggest that the "big blue bump" may 
not be as prominent as previously thought; Laor et al. 1997). The MF87 continuum gives 
rise to line luminosities somewhat larger than observed, as the observed L(Ka) of MCG 
-6-30-15 in the "medium phase" is ?^1.6xl0 41 ergs s _1 (Isawasa et al., 1996). On the other 
hand, the "power-law" continuum has an ionizing luminosity a factor ~ 10 lower than that 
of the MF continuum, for the same optical flux, and it therefore not surprising that fluxes 
are lower than observed. The total Ha luminosity is rather sensitive to the presence of an 
outer halo, but the Fe Ka luminosity is not, as emission in the innermost disk region £ 
20 R g is by far the strongest. Figure 3 shows emissivity profiles for Hi Ha and Fe Ka 
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derived from our models. 

Next we computed the cold Fe Ka line profile following Fanton et al. (1997), and 
using the emissivity law obtained from Model 1. The code developed by Fanton et al. 
computes the trajectories in the full Kerr metric taking into account all the relevant general 
relativistic effects. The computed profile is therefore accurate at any distance beyond the 
event horizon and for disks seen at any inclinations. Integration was carried out in between 
1.25 R g ;$ r ^ 20 R g . Figure 4 shows the results of our fits to the MCG-6-30-15 data 
(Tanaka et al. 1995) where we were able to satisfactorily reproduce the Fe Ka profile. 



3.2.2. Hi Ha 

The vertical stratification of density implies a variety of physical conditions within 
the disk. This has the important implication that the disk is able to produce Balmer 
line emission at almost any distance from the continuum source, in regions (c) and (d) of 
§3.2. 1| . The Ha emissivity law for Model 1 can be approximated with a shallow power law 



(p = -0.71) for 10 2 R g < r < 10 4 R g and with a steeper law (p=-2.33) for r > 10 6 
R g . This emissivity law makes the production of symmetric single-peaked profiles very 
likely, as the emissivity favors emission from the outer regions of the disk, where relativistic 
effects are negligible. e(r) x r 2 peaks at r ~ 10 4 R g , a value that is in agreement with 
the average distance for Balmer line emission deduced from reverberation mapping. Fig. 
5 shows superimposed Ha and disk profiles based on model 1 for 15 Seyfert l's included 
in the Nandra et al. (1997a) sample. Model profiles were computed following Chen & 
Halpern (1989) rather than Fanton et al. (with the emissivity law of Model 1) as their 
weak field approximation is adequate for the relativistic treatment at large radii. Best fits 
were achieved using the downhill simplex algorithm by Nelder & Mead (1978). The outer 
radius of the emitting region was allowed to vary within a narrow range, 10 5 R g £ r out ;$ 
10 6 Rg (the distance at which disk fragmentation due to self gravity is expected to become 
relevant). An additional free parameter is a scale factor. In Table 6 we report the results 
of Hi Ha disk model fitting. The format is: Column 1 and 2 - IAU code and common 
name; Column 3 - r out ; Column 4 - Hi Ha disk inclination from our model 1 fits; Column 5 
- the reduced x 2 ; Column 8 - Column 6 and 7 inclination and reduced x 2 f° r Model 2 fits 
(not shown in Fig. 5); Fe Ka disk inclination from Kerr model fits (Nandra et al. 1997). 
Ha model fits of 3C 120. Statistically good fits are the ones with \ 2 ~ 1- The formal 
uncertainty in the estimate of the inclination is typically a few degrees (it is undefined if 
the fit is not statistically good), but the actual uncertainty is probably much larger. The 
distribution of scattering matter affects the estimate of the inclination, since it change the 
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turning point in the slope of the emissivity laws (Fig. 4). Model 1 and Model 2 should be 
regarded as limiting cases for the emissivity law. 



3.2.3. Main Results and Implications 

Our attempt to constrain the disk models has met with mixed success. We summarize 
here the main achievements and problems. 

• There is poor agreement between the disk inclination derived from X-ray and optical 
lines for MCG-6-30-15. This is an especially disturbing case. The extremely narrow 
Hi Ha profile requires small inclination while the broad Fe Ka profile constrains 

it to pa 30°, even in the case of a Kerr black hole with a~M. Disk inclination is the 
only free parameter for our model because the inner and outer emitting radii are 
constrained by the emissivity law. ~ 7 objects studied by Nandra et al. (1997a,b) 
have derived Ka disk model inclinations well enough constrained for comparison. The 
small number precludes a reliable statistical analysis. Even including some objects 
with uncontrained inclination from Fe Ka data, we find a statistically good fit for 
Hi Ha and inclination consistency between Hi Ha and Fe Ka only for NGC 3227, 
NGC 3783, Mark 841, NGC 4593, NGC 7469 and NGC 5548. 

• Disk models do not reproduce most Balmer line profiles of Seyfert galaxies. This is 
most likely true for all models computed in this paper. There is no known physical 
reason to limit the emitting radius of Hi Ha except for the outer radius at which 
the disk fragments (10 5 - 10 6 R g ). Integrating within this limit, line cores are well 
reproduced only when they are symmetric and unshifted. The shape and asymmetry 
of line wings also does not agree with disk model predictions (Romano et al. 1996). 

• The comparison between observed and disk model Hi Ha profiles suggests a simple 
variant of the old cloud model (see e.g. Mathews & Capriotti 1985, and references 
therein) capable of explaining the new data: a flattened ensemble of dense (n c ~ 10 10 
to ~ 10 13 cm -3 ) clouds, co-axial with the accretion disk, spiralling inward (or outward) 
because of drag forces, on a time-scale much shorter than the radial drift time scale 
for an a-disk. Cloud production might occur in a wind from the disk surface, or the 
wind itself could be a source of Hi Ha and CivA1549 line emission (e. g., Murray & 
Chiang 1997). 

The real challenge for any BLR model is not the ability to reproduce one single 
profile, but rather to account for the distribution of line width, shift and asymmetry 
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observed in any sample of AGN. No model proposed until now has achieved this goal. 
Asymmetries may be produced because the clouds are optically thick to the lines they 
emit (i.e. the lines we observe come mostly from the illuminated face of the cloud), 
and because of the spiralling path. If the system of clouds is observed nearly face on, 
then we have symmetric, unshifted, profiles, and the narrowest of all observed profiles. 
Larger asymmetries are expected if the line is broader (as observed) : in this case line 
emitting clouds have a larger component of their non-circular velocity along the line 
of sight. 

The Fe profile shape may be complicated by a significant and variable contribution 
from hot Fe. This may help account for the large observed EW's and provides an 
additional free parameter for adjustment to the profile shape. In particular it might 
account for the wing on the high energy side of the line observed in sources like 
MCG-5-23-16 (Weaver et al. 1997). 



4. Other Models for Fe Ka Emission 

4.1. Major Challenges for Disk Models 

An important uncertainty, that affects our interpretation of the Fe line origin, is 
whether it arises from a single emitting region. Almost all past disk model fits, including 
the previous section, have treated the Fe Ka line as a single component. There are several 
forms of evidence to support the hypothesis that the line is composite. 

1) Many Seyfert 1 spectra show a narrow Fe Ka peak very close to the rest energy of 
6.4 keV. This is not a problem for a wide range of disk models and is, in fact, consistent 
with the expectation of a Doppler-boosted blue peak from an illuminated disk viewed at 
intermediate inclination. The fact that the blue peak always occurs at 6.4 keV would 
just be a coincidence. However the narrow feature can also be reasonably ascribed to an 
unshifted independent component. This possibility is strengthened because some of the 
best line profiles, including MCG-6-30-15, show a smooth high energy (blue) wing. Nandra 
et al. (1997) present composite spectra (see solid curve in Figure 6) made both with and 
without the two "best" examples, MCG-6-30-15 and NGC 4151. Both composite spectra 
show a similar blue wing on Fe Ka suggesting that it is a common property of Seyfert 1 
line profiles. The range of disk models explored in the literature almost always show a sharp 
drop on the blue side of the profile which is inconsistent with this blue shoulder. The recent 
detection of MCG-5-23-16 (Weaver et al. 1997 - not included in the Nandra et al. 1997 
study), exacerbates this problem showing perhaps the strongest blue shoulder yet observed. 
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An obvious decomposition of the line into two symmetric Gaussian components is: (1) a 
narrow component centered at E n =6.4±0.05 keV with FWHM=0.25 keV and (2) a broad 
component centered at Eb=5.9±0.1 keV with FWHM=1.6 keV. The blue wing appears to 
be a natural extension of the second component under the first rather than an independent 
feature. 

2) Recent variability results for the best studied Seyfert 1 MCG-6-30-15 show that 
the narrow "blue" (6.4keV) and broad red components vary independently (Iwasawa et al. 
1996b). The blue peak responds quickly to continuum changes while the red peak responds 
more slowly or remains relatively constant. This is again consistent with the notion that 
we are observing two independent features. The variability studies add another problem for 
disk models. Pre-ASCA variability studies (Mushotzky et al. 1993) already found potential 
problems with disk models because the Fe Kq line did not show the relatively short time 
scale response to continuum changes that is expected from simple disk models. This issue 
is discussed in a companion paper (Sulentic et al. 1997). 

3) Many objects show breaks in the red wing that could indicate discrete components. 
Only higher s/n observations will decide if the red wing is continuous and consistent, in 
detail, with disk models. The recent Fe Ka detection in MCG-5-23-16 shows clear evidence 
for two or three components (one at 6.4 keV as discussed above) (Weaver et al. 1997). 
Model fits to the broad feature in MCG-5-23-16 leave a significant residual between 5-6 keV 
and possibly also on the blue side. A 5.5 keV component was also discussed for NGC 4051 
(Guainazzi et al. 1996) and 3C 109 (Allen et al. 1997) however the latter reference suggests 
that this may be related to an instrumental artifact. 

4) Finally there is the problem raised by unification models that view (at least some) 
Seyfert 2 galaxies as edge-on Seyfert l's where the BLR is obscured by a dust torus (see 
also Turner et al. 1997 and Nandra et al. 1997). Seyfert 2 galaxies are most probably an 
heterogeneous class, but there is evidence that Seyfert 2 with Fe Ka detections may be 
genuinely obscured Seyfert 1 (see e. g., Dultzin Hacyan 1997, and private communication). 
Emission from an inclined disk would produce broad profile which are easily distinguishable. 
Only Seyfert 2 sources IRAS 18325-5926 and MCG-5-23-16 can be reasonably well fit with 
disk models. The red asymmetry and blue-shifted blue peak are consistent with a reasonably 
high inclination (i=40-50°: Iwasawa et al. 1996a) as expected for a Seyfert 2. The current 
status of disk emission models for Fe Ka lines in Seyfert 2 objects hearkens back to the 
situation for a disk signature in optical broad lines as discussed in the introduction. If 
we see an obvious disk signature in IRAS18325-5926 (and MCG-5-23-16) then where is 
the corresponding signature in other Seyfert 2 objects? Why we do not see an edge-on 
disk profile in NGC 1068 (this object may be optically thick out to 20 keV: Mushotzky et 
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al. 1993), NGC 3147 and 6552 as well as Markarian 3 and other Seyferts 2's observed by 
ASCA? The lack of coherence in the model fits to different AGNs of the same subclass is 
already disturbing. The large equivalent width of the IRAS18325-5926 detection (and now 
also MCG-5-23-16), if from a disk, suggests that many other inclined (but not edge-on) 
disks should be detectable. 



4.2. Alternatives to Fe Ka Disk Emission 

Alternative models for the complex emission profile near 6.4 keV (or just the broad 
component) must avoid the difficulties that were summarized in the preceding section. They 
must also satisfy the following general requirements: (a) produce strong Fe Ka emission 
(EW~100-200eV), (b) broaden it (a ~0.6keV), (c) red shift the bulk of the emission 
and (d) to a lesser extent, symmetrize the line (assuming that the narrow and broad 
components are independent). It must also account for the FWHM difference between 
Fe Ka and the UV/optical HIL/LIL if a BLR unification is attempted. The first task can 
be accomplished with an ensemble of emitting clouds provided that the column density 
is high enough and the covering factor is in the range 55-85% (Sivron and Tsuruta 1993; 
Nandra & George 1994). The strength of the line is a problem for disk and non disk models 
alike. A photoionized medium is one of the solutions previously proposed for this problem 
(e.g. Tanaka et al. 1995). The line profile can be broadened most easily by increasing the 
velocity dispersion of the cloud ensemble but possible constraints imposed by the sharpness 
of absorption edges (Fabian et al. 1995) motivate us to look for solutions that minimize 
Doppler broadening. Several other possibilities have been considered, some of which could 
also introduce the required red displacement that is observed (e.g. Fabian et al. 1995) by 
(1) Comptonization, (2) transverse Doppler effect and (3) gravitational redshift. 

The difference in average profile FWHM between Fe Ka and Hi Ha greatly constrains 
non-disk models where Fe Ka arises in the optical broad line clouds. Our previous work 
on modeling optical line profiles with bi-conical emission structures (as alternatives to disk 
models; Zheng, Binette, & Sulentic 1990) leads us to consider that possibility also for Fe 
emission. A bi-cone model might provide an explanation for the strong changes observed 
in the narrow component of MCG-6-30-15. As the approaching side of the bi-cone should 
be the one closest to the observer, one could see variations following continuum variability 
with essentially zero time delay. This is especially true if the opening angle of the cone 
is small. Dense clouds (n c £ 10 13 cm~ 3 are a likely producer of Fe Ka, and a continuum 
reprocessor as well (Guilbert & Rees 1988; Lightman & White 1988) giving a flattening of 
the X-ray spectrum similar to the situation for a disk reflection component. Such dense 
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clouds should not be a strong source of Balmer line emission. 



4.2.1. A First- Order Model 

We made a first-order test of the possibility that the Fe Ka line emission is produced 
in a cloud ensemble surrounding the central source. Earlier work on a cloud model was 
discussed in the previous section. Our extension of this work is intended to show only 
that a biconical cloud distribution can reproduce the most complex Fe Ka line profiles. 
Options include modeling: (1) the entire Fe line profile for MCG-6-30-15, (2) the entire 
average profile (Nandra et al. 1997) and (3) the broad redshifted part of the average profile. 
The first two options are essentially equivalent so we focus on MCG-6-30-15 as an extreme 
example of option 2. This source has perhaps the best defined and most studied profile of 
any Sy I detection. Options 1 and 2 represent the greater challenges because it is certainly 
easier to reproduce a Gaussian profile with a cloud ensemble model. 

We initially assume that the red component is produced by clouds at rest at a distance 
r moving at speed v from the central source. The gravitational redshift formula: 



A/A' = [1 - (v/c) cos0]/[yj(l - v 2 /c 2 ){l - 2GM/{c 2 r))} 

(where A is the rest frame wavelength, A' is the shifted wavelength, 9 is the angle between 
the line of sight and the cloud velocity, M the mass of the central black hole) then allows 
us to convert the observed profile into an intensity distribution I(r) contributed by shells 
at various radii from the central source. In the case of the redshifted part of the average 
profile the emission peaks at r = 12 R g with a sharp drop towards smaller radii and an 
I(r)oc r^ 1 tail in the opposite direction. The higher redshift of the MCG-6-30-15 profile 
implies even smaller distances from the central source. The bulk redshift of the line is a 
consequence of the fact that the majority of the emission arises very close to the black hole. 
This conclusion is true for any kinematical model, including emission from an accretion 
disk, and does not depend on the specific geometric assumptions. 

Unfortunately it is hard to explain why the clouds at only a few gravitational radii 
from the black hole would be at rest. The gravitational potential is very steep and 
the radiation pressure is unlikely to be important because we are in the sub-Eddington 
regime. A marginally bound particle passes 10 R g with a velocity of 0.45 c during radial 
infall. With motion being dynamically important one must make assumptions about the 
geometry and energetics. Our first-order model considers clouds with predominantly radial 
motion. Angular momentum tends to circularize any cloud distribution therefore making it 
impossible to explain the double peaked line profile observed in MCG-6-30-15 (unless a very 



-20- 



specific planar geometry is assumed as in the case of an accretion disk). We constructed 
models for both radial infall and outflow. As we show below inflow is clearly preferred by 
model fits to the observations. This is consistent with the arguments based on energetics. 
It seems very hard to accelerate the matter outwards in a region only a few R g from the 
black hole. 

The assumptions of the model are as follows: (i) Clouds are in radial bound orbits 
defined by their apastron distance R where the clouds would be at rest. The profile 
is reproduced from a Monte Carlo simulation of individual cloud contributions. The 
distribution of distances r of the clouds is calculated from the time dt clouds are spending 
between shells r and r + dr. (ii) The cloud distribution is limited by the inner radius 
r in (we assume r in = 6 R g ) and the outer radius r out < R (r out = 15 R g was assumed 
throughout), (iii) The clouds are confined to double cones with a half-opening angle i c . 
(iv) The inclination of the line of sight towards the cone axis is i. (v) Emission in the 
comoving frame of the cloud is isotropic in the optically thin case, and confined to the 
inward looking hemisphere for the case of optically thick Fe K« emission (both cases are 
feasible, see Nandra and George 1994). The intensity of emission in the comoving frame 
does not depend on distance from the central source. 

Cone models often suffer from a large number of assumptions involving geometry and 
emissivity laws. Our fitting is only illustrative so we do not pretend to explore the whole 
parameter space here. Instead, by using the model outlined above, we keep the number of 
parameters reasonably low: only the apastron distance and the inclination were fit, with 
the rest being assumed or of second order importance. This is similar to accretion disk 
models where the critical parameters are the inclination and the distance of its inner edge. 

Figure 6 shows the results of our model fits. The thick histogram shows an example 
of the line profile we get for radial inflow of optically thin clouds and the thin histogram 
shows the optically thick counterpart. The inclination (i = 20°) is intermediate and the 
cone is wide (i c = 50°). The value of apastron distance R is important. Here, as in all other 
attempts, we were forced to assume that the velocity of the clouds approaches zero at the 
outer edge of the cloud distribution which is still fairly close to the BH (this is consistent 
with the CLOUDY calculation presented above). Although on one hand a consequence of 
this assumption is that thesse clouds do not reach large distances from the central source 
that would blur the relatively sharp observed O VII edge attributed to the warm absorber 
(Fabian et al. 1994a, 1995), on the other hand it is difficult to explain why the velocity 
would not be higher, i.e. closer to that of marginally bound infalling clouds. The accretion 
disk scenario escapes extreme velocities by invoking nearly face-on orientation. In a bi-cone 
model such a projection effect seems unlikely unless an inclination nearly perpendicular to 
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the cone axis is assumed. In this case the contributions from both cones merge in a single 
Gaussian-like redshifted component. It is easy to fit the redshifted Gaussian component 
of the average profile this way, but not the well defined peak and red shoulder of the 
MCG-6-30-15 profile. We do not attempt to explain the low radial velocity of clouds at 
~ 15R g here. Possible explanations include a contribution from an ion torus or from a 
geometrically thick inner edge of a surrounding accretion disk. Mathews (1993) was able to 
"bounce" clouds around an equilibrium radius (where gravity was balanced by radiation 
forces) about an order of magnitude larger than the one considered here. 

It would be satisfying to fit the average profile using the same parameters as used 
for the MCG-6-30-15 fit, varying only the inclination to the line of sight. Unfortunately 
we can not envision a model where this is true, accretion disk models included. In the 
case of a bi-cone the higher inclination angle causes the peaks to merge so that the narrow 
component at the rest energy must have a separate origin (an interpretation that we favor). 
A convenient way to fit both components of the average profile with a single zone model 
would involve the assumption of varying optical thickness of the emitting clouds. The thin 
lined histogram in Fig 5 uses the same parameters as the MCG-6-30-15 fit, but assuming 
the optically thick regime. The low energy peak is weaker as a consequence of the fact that 
the optically thick inward moving clouds on our side of the BH emit most of light away 
from us, while we still see the bright faces of the clouds that lie on the other side and move 
toward us. The observed average profile lies between both histogram curves and therefore 
it is reasonable to conclude that the full range of observed Fe Ka lines can be explained 
by different optical thickness of the clouds (Nandra and George 1994). 

The model outlined above allows us to distinguish between two kinds of radial motion, 
i.e. inflow and outflow. Some of the photons originating in the far-side cone are intercepted 
by the BH, so the contribution from the far-side cone is weaker. The blue part of the profile 
gets weaker if inflow is considered, while the red part is weaker in the case of outflow. In 
all of our Monte Carlo simulations it was hard to avoid too much flux in the high energy 
wing of the line. The contribution from clouds close to the black hole moving towards the 
observer was acceptably low only for the case of inflow (though the high energy flux is still 
higher than observed, cf. the thick histogram in Fig. 5). The distinction between both 
directions of motion is even more obvious in the explanation of the difference between the 
MCG-6-30-15 and the average profile as outlined above. If we deal with outflow it is the 
high energy peak that gets weaker when we switch from the optically thin to the thick 
regime, i.e. opposite to the case observed. 

We conclude that a simple double cone model has the potential to explain both 
the MCG-6-30-15 and the average Fe Ka profile allowing only for a variation of optical 
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thickness of the emitting clouds. This can be further used to explain the observed variations 
in some objects. This is encouraging, still a physical explanation of the origin of low velocity 
clouds at ~ 15R g is clearly needed but is beyond the scope of this paper. 



5. Summary of Results & Implications for AGN Models 

The main results of this investigation are as follows: 

1. The bulk of Fe emission near 6.4 keV cannot be directly associated with the optical 
BLR. Fe Ka and Hi Ha lines show significantly different width and are therefore 
not emitted in the same BLR region in radio-quiet objects, at least not at the same 
distance from the central continuum source. 

2. A simple illumination model suggests that both Fe Ka and Hi Ha could be emitted 
by an accretion disk, although the bulk emission for the two lines would arise at 
widely different radii (Fe Ka < 100 R g ; Hi Ha~ 10 3 R g ). 

3. While our illuminated disk model is able to predict Hi Ha profile widths and emitting 
radii consistent with observational results, the observed and model profile shapes 
show poor agreement. The Fe Ka emission line data are also inconsistent with a 
single emissivity law. There is poor agreement between disk inclination values derived 
for Fe Ka and Hi Ha. Observed Hi Ha asymmetries and shifts demand a larger 
radial drift velocity than compatible with Keplerian disk models. 

4. Evidence is accumulating to suggest that the Fe Ka line may be composed of (at 
least) two independent components. If the broad component is independent and 
also shows Gaussian symmetry as suggested by the Nandra et al. (1997b) composite 
spectrum, there are no obvious disk models consistent with it. Physically, there are 
viable alternatives. The narrow component could be due to BLR emission with a 
smaller contribution from the warm absorber, while the broad component could be 
redshifted coronal emission. 

5. First-order models involving a biconical distribution of emitting clouds are able to 
reproduce the observed Fe Ka spectra for AGN. They can also easily account for the 
observations that raise serious problems for the disk emission scenario. 
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Fig. 1. — Upper: Spectrum of Arpl02B in the region of Hi Ha. Taken from a 35 minute 
exposure at the Asiago Observatory 1.82 m. with a dispersion of 120 A/mm (March 26, 1989 
UT=0254). Lower: Same spectrum degraded to the sampling rate of an ASCA observation 
of Fe Ka (after narrow lines subtraction). Solid line: disk model profile fitted to the Ha 
profile by Chen & Halpern (1989); dot-dashed line gaussian profile. 



Fig. 2.— FHWM of Fe Ka versus FWHM of Hi Ha, in km s" 1 . Filled circles: SIS ASCA 
data; filled squares: BLRG; circled filled circles: best candidate disk emitters; the star is 
NGC 1068 (broad Ha component visible in polarized light). The dot-dashed line is the line 
of equal Fe Ka and Hi Ha FWHM. 



Fig. 3. — Plots of Hi Haand Fe Ka surface emissivity, in units of ergs s _1 cm -2 , derived 
from our illumination models as a function of disk radius. Full line: "cold" Fe Ka (iron 
less than 16 times ionized) emissivity; dashed line: "hot" Fe Ka emissivity; dotted line: 
Hi Ha emissivity. See text and Table 4 for model parameters. Labels identify the model as 
in Table 4. "Hot" iron emission has been plotted for model (2) only for clarity. Ha surface 
emissivity has been plotted for r £ 100 R g . 



Fig. 4. — Fe Ka line profile for MCG-6-30-15 as obtained by Tanaka et al. 1995. The solid 
line is our best Kerr disk profile model fit with r; n m 1.25R g , r out ~ 20R g , i = 30° and 
power-law emissivity with spectral index p = —1.8. 



Fig. 5.— Optical Ha line profiles of 15 AGN considered in this paper, after subtraction 
of the underlying continuum. Vertical scale is specific flux (ergs s" 1 cm" 2 A -1 x 10 15 ); 
horizontal scale is wavelength in A. The thick lines show model profiles computed following 
Chen & Halpern (1989). Dotted lines indicate the pure broad line (narrow lines substracted) . 
See Table 6 for model parameters. 



Fig. 6. — Observations and theoretical models of the Fe Ka line. Points: observed data for 
MCG-6-30-15, dashed curve: accretion disk fit (Tanaka et al. 1995), solid line: observed 
averaged profile for 16 sources (Nandra et al. 1997). (a) thick histogram: double-cone 
optically thin radial inflow of 10 4 clouds (i = 20°, i c = 50°, r = 16R g : see text), thin 
histogram: same for optically thick inflow. 
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Table of ASCA FeK line detections 
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NOTE. — Col. (1) - abbreviated (2000) position taken from the 7th edition of the Veron & 
Vcron catalog; Col. (2) - a common name for the source as used in the ASCA-related literature; 
Col. (3) - AGN type; Col. (4) - optically determined redshift; Col. (5), (6), and (7) - the derived 
Fe Kq line peak centroid energy, profile Gaussian width (kcV) and EW (cV), respectively; Col. 
(8) - refcrcnce(s) to the source of the Fe Kq parameters given in the three previous columns. 
*: objects with clear multicomponcnt FeK line profiles, see references for details. 

REFERENCES. — 1: Tanaka ct al., 1995; 2: Yaqoob ct al., 1995; 3: Mushotzky et al., 1995; 
4: Nandra et al., 1997; 5: Eraclcous et al., 1996; 6: Iwasawa et al., 1996a; 7: Weaver et al., 
1995; 8: Ishisaki et al., 1996; 9: Turner et al., 1996; 10: Mihara ct al., 1994; 11: Pounds ct al., 
1994; 12: Nandra et al., 1997; 13: Ptak ct al., 1994; 14: Guainazzi et al., 1996; 15: George et 
al., 1995; 16: Weaver et al., 1996; 17: Pounds et al., 1995; 18: Reynolds et al., 1995; 19: Ptak 
et al., 1996; 20: Iwasawa et al., 1994; 21: Ueno ct al., 1994; 22: Hayashi et al., 1996; 23: Allen 
et al., 1997; 24: Yaqoob et al., 1994a; 25: Fukuzawa et al., 1994; 26: Yaqoob et al., 1994b; 27: 
Makishima ct al., 1994; 28: Weaver et al., 1997; 29: Fabian et al., 1994; 30: Done et al., 1996, 
ApJ; 31: Elvis et al., 1994. 
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B&C 


600 gr/mm 


4 


J1210+394 


NGC 4151 


Mar. 24, 


1989 


00:42 


1800 


Asiago 


1.82m 


B&C 


600 gr/mm 


4 


J1239-053 


NGC 4593 


Apr. 2, 


1990 


06:03 


3600 


ESO 


1.52m 


B&C 


1200 gr/mm 


2 


J1335-342 


MCG -06-30-15 


Feb. 3, 


1995 


11:39 


2400 


SPM 


2.12m 


B&C 


1200 gr/mm 


2 


J1504+104 


Mark 841 


May 21, 


1993 


04:13 


3300 


ESO 


1.52m 


B&C 


1200 gr/mm 


2 


J 1942 -103 


NGC 6814 


May 23, 


1993 


06:37 


1800 


ESO 


1.52m 


B&C 


600 gr/mm 


4 


J2044-107 


Mark 509 


May, 25, 


1993 


09:33 


3600 


ESO 


1.52m 


B&C 


600 gr/mm 


4 


J2303+088 


NGC 7469 


Aug. 20, 


1990 


10:28 


1200 


SPM 


2.12m 


B&C 


600 gr/mm 


4 


J2304-086 


Mark 926 


Sept. 26, 


1989 


03:12 


3000 


ESO 


1.52m 


B&C 


1200 gr/mm 


2 



NOTE. — Col. (1) - abbreviated IAU designation; Col. (2) 
of observation; Col. (5) - exposure time in seconds; Col. (6) 
grating; Col. (10) - resolution FWHM (in A). 



other name; Col. (3) - 

observatory; Col. (7) 



date of observation; 
telescope; Col. (8) - 



Col. (4) - universal time 
spectrograph; Col. (9) - 



TABLE 3 
Broad Ha Line Profiles 



IAU Code 

(1) 


Name 
(2) 


S. T. 
(3) 


EW 

(4) 


FWHM 

(5) 


A(FWHM) 
(6) 


Vr(l/2) 

(7) 


Av r 

(8) 


Ref. 

(9) 


J0006+202* 


Mark 335 


Syl 


490 


1600 


300 


300 


100 


1,2 


J0123-588* + 

J0228+313 

J0242-000* 


Fairall 9 
Mark 1040 
NGC 1068 


Syl 

Syl 
Sy2 


420 
305 


4550 
2200 
4100 


500 
300 
1000 


410 



200 
200 


3 

4,5 
6 


J0433+053* + 
.10558-383 
J0913+409 
J0955+690 


3C 120 

EXO 0556.3-3820 
IRAS 0910 
NGC 3031 


EG 

Syl 

QSO 

Syl 


470 
760 

26 


2400 
4250 

3400 


300 
400 

350 


110 
-160 

-140 


100 
100 

200 


7 

8 

9.10 


,11023+198* + 
,11034+396 


NGC 3227 
ZW 212.025 


Syl 
Syl 


285 
200 


4000 
1800 


400 
200 


70 



100 
200 


7 
11 


,11106+725* + 


NGC 3516 


Syl 


150 


5300 


500 


-270 


100 


12 


,11119+213+ 
,11139-377* + 
J1203+445* 
J1210+394* 


PG 1116+215 
NGC 3783 
NGC 4051 
NGC 4151 


QSO 

Syl 

Syl 

Syl. 5 


340 
440 

240 


2750 
3050 
1000 
6400 


400 
300 
200 
600 


60 

40 
100 
-350 


100 
100 
100 
100 


7 

3 
5 

12 


J1214+140+ 
J1218+298* 
J1229+020 


PG 1211 + 143 
Mark 766 
3C 273 


QSO 

Syl 

QSO 


390 
*195 
613 


1700 
1600 
3400 


300 
200 
400 


-30 


100 


100 
200 
100 


7 

13,14 

15,2 


,11239-053* + 


NGC 4593 


Syl 


250 


3400 


400 


-180 


70 


12 


,11335-342* + 


MCG -6-30-15 


Syl 


160 


2200 


200 


50 


100 


12 


,11349-303* + 


IC 4329A 


Syl 


515 


5000 


400 


330 


150 


16 


J1417+251* + 


NGC 5548 


Syl 


540 


5450 


300 


-590 


200 


3 


J1504+104* + 
J1842+797 


Mark 841 
3C 390.3 


Syl 
KG 


470 
450 


5000 
12000 


400 
1000 


-220 
300 


100 
150 


12 

12,7 


,11942-103* + 


NGC 6814 


Syl 


50 


5600 


1000 


360 


200 


12 


J2044-107* + 


Mark 509 


Syl 


550 


3200 


400 


-200 


100 


12 


J2303+088* + 
J2304-086* + 


NGC 7469 
Mark 926 


Syl 
Syl 


380 
800 


2950 
7650 


300 
1500 


320 
790 


100 
300 


12 
12 



NOTE. — Col. (1) - IAU code name; Col. (2) - source name; Col. (3) - Seyfert type; Col. (4) - EW in A; 
Col. (5) - FWHM in km/s; Col. (6) - uncertainty at 2xcr confidence level; Col. (7) - centroid line shift (at 
half maximum) in km s~ ; Col. (8) - line centroid uncertainty; Col. (9) - reference code. *: objects studied by 
Nandra et al. 1997; f : Ha profile shown in Fig. 2; : broad component revealed in polarized light; : estimated 
as EW(H/3) X Balmcr decrement. 

REFERENCES. — 1: Shudcr, 1981; 2: Cassidy & Rainc, 1993; 3: Marziani ct al., 1996; 4: Rafanelli & Schulz, 
1983; 5: Ostcrbrock & Shudcr, 1982; 6: Antonucci & Miller, 1985; 7: G. Stirpe, private communication & 
unpublished thesis work; 8: Morris & Ward, 1989; 9: Filippenko & Sargent, 1988; 10: Ho, Filippenko & Sargent, 
1995; 11: Puchnarewicz ct al., 1995 12: this work: 13: Ostcrbrock & Poggc, 1985; 14: Gonzalez-Delgado & 
Perez, 1996; 15: Morris & Ward, 1988; 16: Marziani, Calvani, & Sulentic, 1992. 



TABLE 4 
Explorative Photoionization Computations 



Source 



log L(Fe Ka) 



Fe Ionization 



log 10 (Hl Ha) 



(1) 


A(Fe) q 

(2) 


lOx A(Fe) q 
(3) 


(4) 


(5) 


BLR 

Warm Absorber 


39.65 
39.22 


40.50 
40.24 


II 

IX 


41.5 
40.0 


Coronal gas 


£ 35.7 


Z 38.2 


XXVII 


32.4 


Dusty torus 


5 38.85 


£ 39.28 


IV 





NOTE. — Col. (1) - possible origin for Fe Ka line emission; Col. (2) and (3) - the 
Fe Ka luminosity computed in the case of gas with solar chemical abundances and 
of a factor 10 iron overabundance; Col. (4) - the dominant ionization stage of iron; 
Col. (5) - the expected Hi Ha luminosity. 

a L(Hl Ha) cannot be reliably predicted by an optically thin torus model. 



TABLE 5 
Disk Model Results for MCG -06-30-15 



Model r ha lo,i 
(1) (2) 



rwo.o. 
(3) 



(4) 



s Halo 
(5) 



Continuum 

(6) 



5io L(Fe Ka) 

(7) 



Feiih. 



,t/FEII 00 

(8) 



log 10 L(Hi Ha) 

(9) 



(1) 


1. 


100. 


0.25 


(2) 


1. 


100. 


0.25 


(3) 


1. 


100. 


0.25 



0.25 



0.25 



MF97 


41.86 


Z 0.3 


41.62 


MF97 


41.36 


0.1 


40.63 


•owcr-law 


40.56 


negligible 


39.37 



NOTE. — Col. (1) - identification number for each model; Col. (2) (3), (4) and (5) - the inner and outer radii, the electron 
scattering optical depth for the inner sphere and the outer halo; Col. (6) - the shape of the continuum; Col. (7), (8) and (9) - 
the logarithms of the Fe Ka luminosity from cold iron emission, the ratio L(Fe Ka)Hot/L(FE Ka)cold, and the log(L(Hl Ha)) 
respectively. 



TABLE 6 
Broad Ha Disk Model Fitting Results 



IAU Code 

(1) 


Name 

(2) 


Rout 

(3) 


io P t[l] 
(4) 


x 2 [i] 

(5) 


io P t[2] 
(6) 


X 2 [2] 
(7) 


i F <. 
(8) 


J0123-588 


Fairall 9 


l.OxlO 5 


36° 


6.16 


15° 


0.90 


89° 


J0433+053 


3C 120 


1.9x10 s 


34° 


1.00 


9° 


1.58 


88° 


J1023+198 


NGC 3227 


1.0x10 s 


36° 


4.88 


12° 


0.81 


21° 


J1106+725 
J1139-377 


NGC 3516 
NGC 3783 


1.5x10 s 
2.0x10 s 


36° 

50° 


9.96 
1.00 


20° 
12° 


1.63 
2.25 


26° 
26° 


.11210+394 


NGC 4151 


1.5x10 s 


50° 


9.44 


25° 


1.38 


24° 


J1239-053 


NGC 4593 


1.0x10 s 


38° 


1.36 


14° 


0.76 


0° 


J1335-342 


MCG-6-30-15 


l.OxlO 6 


13° 


0.46 


9° 


2.95 


34° 


J1349-303 


IC 4329A 


1.0x10 s 


36° 


13.6 


12° 


5.24 


10° 


J1417+251 


NGC 5548 


1.0x10 s 


37° 


7.64 


18° 


0.25 


10° 


J1504+104 
J1942-103 


Mark 841 
NGC 6814 


1.0x10 s 
1.0x10 s 


38° 
38° 


4.36 
0.67 


17° 
19° 


0.83 

0.1 


26° 


J2044-107 


Mark 509 


1.0x10 s 


35° 


4.64 


12° 


7.08 


89° 


.12303+088 


NGC 7469 


9.0x10 s 


50° 


1.04 


10° 


0.97 


20° 


J2304-086 


Mark 926 


1.0x10 s 


38° 


35.0 


26° 


2.77 


26° 



NOTE. — All objects studied by Nandra et al. 1997 and with Ha profile shown in Fig. 
5. Col. (1) and (2) - IAU code and common name; Col. (3) - disk outer radius; Col. 
(4) - Hi Ha disk inclination from our model 1 fits; Col. (5) - the reduced \ j Col. (6) - 
Hi Ha disk inclination from our model 2 fits; Col. (7) - the reduced x ! Col. (8) - Fe Ka 
disk inclination; data taken from Nandra et al. (1997) for Kerr geometry with reflection: 
in cases where multiple observations were presented, we chose the inclination with the 
smallest uncertainty. 



